Electronic structure of half-metallic double perovskites 
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We present the self-interaction corrected local spin density (SIC-LSD) electronic structure and 
total energy calculations, leading also to valencies of the ground state configurations, for the half- 
metallic double perovskites such as Sr2FeMo06, Ba2FeMo06, Ca2FeMo06, and Ca2FeRe06. We 
conclude that the Fe and Mo (or Re) spin magnetic moments are anti-parallel aligned, and the 
magnitude of the hybridization induced moment on Mo does not vary much between the different 
compounds. The hybridization spin magnetic moment on Re is of the order of -1.1 fiB, while that 
on Mo is about -0.4 ^b, independently of the alkaline earth element. Also the electronic structure 
of all the compounds studied is very similar, with a well defined gap in the majority spin component 
and metallic density of states for the minority spin component, with highly hybridized Fe, Mo (or 
Re), and oxygen bands. 

PACS numbers: 75.50.Gg, 75.47.Pq 



I. INTRODUCTION 

The interest in double perovskites A2B'B"06, with A 
being alkaline earth (AE) ion, and the transition metal 
sites (B-sites) alternately occupied by different cations B' 
and B" , has really taken off with the paper by Kobayashi 
et al.^, demonstrating a convincing half-metallic be- 
haviour in Sr2FcMo06. The calculated electronic struc- 
ture of this compound shows a gap in the majority spin 
component, while for the minority spin component one 
observes strongly hybridized Fe 3c?(t2g), Mo Ad{t2g), and 
O 2p states at the Fermi level. Their study has also estab- 
lished that the latter compound exhibits tunnelling mag- 
netoresistance (TMR) at low applied magnetic field even 
at room temperature. The double perovskites were dis- 
covered in 1960 0^1^1^1^1^ and seem to be metallic and ferro- 
magnetic for B'=Fe, and B"=Mo or Re. Unlike the man- 
ganites, they are stoichiometric and have large Curie tem- 
peratures, Tc, which makes them very promising for tech- 
nological applications. The Sr2FeMo06 compound has 
Tc of about 450 K, but also BazFeMoOe and CasFeMoOe 
have been shown to have interesting TMR characteristics 
and high T^'s.^ Although SrzFeReOe (see Ref. |^ shows 
similar TMR properties to those of Sr2FeMo06, recent 
studies claim that the Ca-compound, Ca2FeRe06, un- 
like Ca2FeMo06, may be insulatingiS It has been sug- 
gested that the magnetic temperature Tc, and whether 
the ground state is metallic or insulating, vary as A is 
changed from Ba to Sr to Ga,^ Also mis-site disorder 
may have pronounced effect on the magnetic properties 
of these compounds >ii*i^ 

Recent experimentsi^ have established in Sr2FeMo06 
an electronic structure and configuration of five localized 
majority electrons (high spin state) on Fc and one dclo- 
calized electron shared between Mo and the other sites. 
We aim to validate this electronic structure by a more 
systematic study involving two other alkaline earths, Ca 
and Ba, as well as Sr. The aspect of a delocalized Mo 



electron is investigated by replacing Mo by Re and in- 
vestigating the amount of delocalized electrons shared 
between Re and the other sites. 

In this paper we present an application of the self- 
interaction corrected local spin density (SIC-LSD) ap- 
proximation to the double perovskites Ba2FeMo06, 
Ca2FeMo06, SraFeMoOe, and Ca2FeRe06. Specifically, 
we concentrate on the electronic and magnetic properties 
of these compounds and in particular the size and rela- 
tive orientation of the spin moments of Fe and Mo (or 
Re). In all the compounds the SIC-LSD calculations find 
a spin moment of about -0.4 fj,B on the Mo and about 
-1.1 fiB on Re, with an opposite orientation to the re- 
spective spin moments on Fe sites. In addition, we find 
these double perovskites to be half-metallic, with a well 
defined gap in the spin-up density of states, and strong 
hybridization at the Fermi energy between the spin-down 
Fe 3d, Mo M (or Re 5d), and O 2p states. 

The earlier first-principles calculations performed for 
double perovskites had to be done with GGA and im- 
plement theoretical optimisation of oxygen positions to 
obtain the half-metallic groundstateii*i^ Other applica- 
tions were done with the LDA-I-U,-'^^ which similarly to 
SIC-LSD provides better description of localized states in 
solids. The present study is the first SIC-LSD applica- 
tion to double perovskites. Its advantage is that it treats 
both localized and delocalized states on equal footing, 
and therefore should be well suited for describing the lo- 
calized nature of Fe 3d electrons as well as the Fe valency. 

The remainder of the paper is organized as follows. In 
the next section, we briefly present the theoretical back- 
ground of the SIC-LSD approach. In section III, the tech- 
nical and computational details regarding the application 
to double perovskites are elaborated upon. In section 
IV, we present total energy and density of states calcu- 
lations. There we also compare magnetic moments and 
determine valency of Fe in the studied compounds. The 
paper is concluded in section V. 
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II. THEORY 

The basis of the SIC-LSD formahsm is a self- 
interaction free total energy functional, E^^'-'' , obtained 
by subtracting from the LSD total energy functional, 
E^^^ , a spurious self- interaction of each occupied elec- 
tron state ij^dT^i namely 

occ. 

OL 

Here a numbers the occupied states and the self- 
interaction correction for the state a is 

= (2) 

with U[na\ being the Hartree energy and E^^^[na\ 
the LSD exchange-correlation energy for the correspond- 
ing charge density ria and spin density tIq,. It is the 
LSD approximation to the exact exchange-correlation 
energy functional which gives rise to the spurious self- 
interaction. The exact exchange-correlation energy E^c 
has the property that for any single electron spin density, 
fiQ, it cancels exactly the Hartree energy, namely 

C/[nJ+£;,,[fv] -0. (3) 

In the LSD approximation this cancellation does not 
take place, and for truly localized states the correc- 
tion (Eq. (PI) can be substantially different from zero. 
For extended states in periodic solids the self-interaction 
vanishes. Consequently, the SIC-LSD approach can be 
viewed as a genuine extension of LSD in the sense that 
the self-interaction correction is only finite for spatially 
localized states, while for Bloch-like single-particle states 
E^"^ is equal to E^^'^ . Thus, the LSD minimum is also 
a local minimum of E^'^'^' . 

It follows from minimization of Eq. that within 
the SIC-LSD approach such localized electrons move in a 
different potential than the delocalized valence electrons 
which respond to the effective LSD potential. For ex- 
ample, in the case of the double perovskites, five (Fe^+) 
or six (Fe^"'') Fe d electrons move in the SIC potential, 
while all other electrons feel only the effective LSD po- 
tential. Thus, by including an explicit energy contribu- 
tion for an electron to localize, the ab-initio SIC-LSD 
describes both localized and delocalized electrons on an 
equal footing, leading to a greatly improved description 
of static Coulomb correlation effects over the LSD ap- 
proximation. Assuming various atomic configurations, 
consisting of different numbers of localized states, one 
can explore the corresponding local minima of the SIC- 
LSD energy functional of Eq. (^), and determine the 
lowest energy solution and valency. The advantage of 
the SIC-LSD formalism is that for such systems as tran- 
sition metal oxides or rare earth compounds the lowest 
energy solution will describe the situation where some 
single-electron states may not be of Bloch-like form. Re- 
garding the valency, it is determined here as the integer 



number of electrons available for band formation, i.e., 

Nval ^ Z — Ncore ~ NsiC, 

where Z is the atomic number (26 for Fe), Ncore is the 
number of core (and semi-core) electrons (18 for Fe), and 
Nsic is the number of localized, i.e., self-interaction cor- 
rected, states (either five or six, respectively for Fe'^+ and 
Fe2+). 

In the present work the SIC-LSD approach has 
been implementedii within the linear muffin-tin-orbital 
(LMTO) band structure method^^ in the tight-binding 
representationiS, where the electron wave functions are 
expanded in terms of the screened muffin-tin orbitals, and 
the minimization of E^^^ becomes a non-linear problem 
in the expansion coefficients. The atomic spheres approx- 
imation (ASA) has been used, according to which the 
polyhedral Wigner Seitz cell is approximated by slightly 
overlapping atom centered spheres, with a total volume 
equal to the actual crystal volume. 

III. CALCULATIONAL DETAILS 

All the compounds studied here crystallize in ordered 
double perovskite structures. The simplest system is 
Ba2FeMo06 which occurs in the cubic structure, with 
the Fm3m space group. The Sr2FeMo06 compound 
has the body-centered tetragonal bet structure, with 
the I4/mmm space group. Replacing Sr by a smaller 
still Ca element results in a monoclinic structure for 
the Ca2FeMo06 and Ca2FeRe06 compounds, with the 
P2i/n space group. The experimental lattice parameters 
used in the calculations for the first three compounds 
have been taken from Refs. and for the last com- 

pound from Ref. ,9. For the Ca compounds two formula 
units (20 atoms) have been necessary. In addition, for a 
better space filling also eight empty spheres have been in- 
troduced in Re-compound. Concerning the linear muffin- 
tin basis functions, we have used 4s, 4p, and 'id partial 
waves on all Fe atoms, and treated them as low- waves, on 
the oxygens (2s and 2p) and Ca (4s and 3p) in addition 
to the s and p low partial waves, 3ci— partial waves have 
been treated as intermediate. On the empty spheres the 
Is— waves have been considered as low- and the 2p's as 
intermediate- waves. On the Ba (6s, bp, and 5d), Sr (5s, 
4p, and id), Mo (5s, 5p, and Ad) and Re (6s, 6p, and 
5d) atoms, the s— , p—, and d— waves have been treated 
as low and the 4/— waves as intermediate. In all the 
cases, the SIC-LSD calculations have been started using 
the self-consistent LSD charge densities with the Perdew 
and Zunger exchange-correlation potential. Mostly the 
Fe 3d electrons have been treated as localized states (in 
some test calculations we have also treated Mo d's as 
localized), meaning that they have been moving in the 
self-interaction corrected LSD potentials, while all other 
electrons have seen the effective LSD potential. In all the 
cases studied, we have allowed for a number of different 
configurations of SIC 3d states to find the ground state 



3 



solution. Namely, treating five 3d electrons on the Fe's 
as localized, moving in the SIC-LSD effective potentials, 
meant that the Fe ions were trivalent. Similarly, for di- 
valent Fe's, six 3d electrons on the Fe's were considered 
to be localized. By comparing total energies correspond- 
ing to different localized Fe 3d configurations we have 
been able to find the ground state configuration which 
for all the compounds has been that of Fe^+. However, 
other configurations have been energetically close, with 
the Fe^+ being insulating in all Mo-based compounds. 
Although, unlike in the work by Kobayashi et al.ji op- 
timizing oxygen positions in Sr2FeMo06 was not needed 
for obtaining half-metallic ground state, however in or- 
der to study its effect on the magnetic properties and to 
make contact with Ref. for this compound we have 
also performed calculations for the scenario where the 
theoretically optimized positions of Ref. Q have been 
implemented 



IV. RESULTS AND DISCUSSION 

The SIC-LSD spin-polarized total densities of states 
(DOS) per formula unit are presented in Fig. 1 for 
CaaFeReOg, CaaFeMoOg, SrzFeMoOg, and BaaFeMoOg 
compounds. There we also show the O 2p spin polar- 
ized densities of states. The spin polarized densities of 
states of Fe 3d and Mo Ad and Re 5d are given in Fig. 
2, for all the different compounds. First thing to note 
is the similarity between the densities of states, indepen- 
dently of which alkaline earth compound is considered. 
For all compounds we see convincing half-metallic be- 
haviour, with well defined gap at the Fermi energy in 
the majority spin channel, and strongly hybridized Fe 
3d, Mo id (or Re 5d) and oxygen 2p states, in the other 
spin channel. The half-metallic character is reflected in 
the total spin magnetic moments that are all integer, as 
can be seen in Table I, where we summarize the total 
and species decomposed spin moments for all the com- 
pounds together with the respective volumes per formula 
unit. As can be seen in Fig. 1, the valence band consists 
predominantly of O 2p states, but contains a small ad- 
mixture of Mo, Re- and Fe d— states. The unoccupied 
majority Re d bands lie closer to the Fermi energy in 
comparison with the unoccupied majority Mo d states 
in Ca2FeMo06. We can also see that the unoccupied 
Mo d states move closer to the Fermi energy as one tra- 
verses the series in the sequence Ca, Sr and Ba. How- 
ever, the majority spin channel maintains a band gap 
for all these compounds. In the minority spin channel 
these systems are metallic with O 2p, Mo 4d and Fe 3d 
states straddling the Fermi energy. These are the states 
which reduce the Fe magnetic moment from essentially 
5 fiB to 3.87 (CaaFeReOe), 3.76 ^j.b (CaaFeMoOg), 
3.71 fiB (SraFeMoOe) and 3.81 hb (BaaFeMoOg) (see 
Table I). From these DOS one can also expect the oxy- 
gen to acquire a moment. This happens according to 
Table I and it is parallel to the Fe spin moment and 



TABLE I: Total and species decomposed spin magnetic mo- 
ments (in /ifl) as calculated self-consistently within SIC-LSD 
band structure method. In the last row the respective vol- 
umes per formula unit are quoted (in (atomic unit)''). 

Moment Ca2FeRe06 Ca2FeMo06 Sr2FeMo06 Ba2FeMo06 





3.00 


4.00 


4.00 (4.00) 


4.00 




3.87 


3.76 


3.71 (3.65) 


3.81 


Mmo 




-0.40 


-0.43 (-0.35) 


-0.41 




-1.12 








Mae 


-0.02 


0.01 


0.02 (0.02) 


0.02 


Moi 


0.02 


0.10 


0.11 (0.11) 


0.09 


Mo2 


0.01 


0.11 


0.11 (0.11) 




Mo3 


0.11 


0.11 






Volume 


773.8 


777.5 


830.2 


884.0 



is non-neglegible at around 0.1 ^b- However it is the 
Mo and Re sites which are of interest. Substantial spin 
magnetic moments of -0.4 /i^ and -1.1 fiB are found on 
the Mo and Re sites respectively. This results from the 
strong hybridization with Fe in the minority spin chan- 
nel. This also explains the antiparallel alignment of these 
spin magnetic moments with respect to Fe. 

The densities of states in the majority spin channels for 
these four compounds are insulating, since a fully occu- 
pied O p band and five localized Fe d states accomodate 
23 (6x3-1-5=23) valence electrons. Of the remaining 19 
or 20 valence electrons for the Mo and Re compounds re- 
spectively, 18 fill up the minority O p band. This leaves 
one or two electrons for the Mo and Re compounds to oc- 
cupy the minority states which straddle the Fermi level, 
and have Fe d, Mo or Re d and oxygen p character. We 
have also considered a scenario of localising through the 
self-interaction correction one Fe (Mo) d state in the mi- 
nority spin channel describing Fe^-|- (Mo^-|-) respectively. 
In this case the states straddling the Fermi energy which 
accomodate this one electron are pulled down below the 
bottom of the O p valence bands and we also obtain an 
insulating state in the minority spin channel. However, 
these scenarios are energetically unfavourable by 70 mRy 
and 103 mRy for Mo^+ and Fe^"'' respectively. Obviously 
the gain in localisation energy does not compensate the 
loss in hybridization energy. 

It is interesting to see how little the magnetic proper- 
ties of these compounds are affected by the size of the 
alkaline earth atom. There is hardly any variation in the 
size of the Fe and Mo or Re moments with the substantial 
change in volume of the compounds (see Table I). Also 
the reduction in the size of the Fe and Mo spin moments 
for Sr2FeMo06, as a result of implementing theoretical 
optimization of oxygen positions, as given by Kobayashi 
et al.,^ does not seem too important although brings the 
size of our calculated moments close to the values of Ref. 
^ The reduced values are given in parenthesis in the 
fourth column of Table I. Concerning spin moments of 
the other species, the small changes cannot be resolved 
when rounded off to two digits after the decimal point. 
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Finally, we would like to comment on the valency of 
Fe and the ground state configuration for all the com- 
pounds. Regarding energetics, in all SIC-LSD calcula- 
tions the Fe'^+ configuration has been most energetically 
favourable, followed by Fe^+ (energetically unfavourable 
by 53, 61 and 83 mRy in Ba-, Sr- and Ca-compounds, 
respectively), and Fe^+ solution, at least twice as un- 
favourable as the former configuration (specifically, by 
103, 122 and 250 mRy). Moreover, treating also Mo d 
states as localized has been energetically unfavourable, 
and the more states have been localized the more un- 
favourable solutions have been obtained. 



V. CONCLUSIONS 

In summary, the SIC-LSD approach has been applied 
to four different double pcrovskites, and in all cases the 
half-metallic ground state solutions have been obtained. 
This seems to be a rather generic result for double pcr- 
ovskites involving Fe with a filled majority d shell and 
transition metal ions which donate few, one (Mo) to two 
(Re), d electrons. This type of materials could there- 
fore be a fruitful playground for the discovery of new 



spintronics materials. In particular, also the Ca2FeRe06 
has been found to be half-metallic, and not insulating as 
indicated in Ref. 0. The electronic structure for all stud- 
ied compounds appears to be very similar, with a gap 
in the majority spin component, and strongly hybridized 
bands at the Fermi level for the minority spin compo- 
nent. The spin magnetic moments induced by hybridiza- 
tion on the Mo sites have been of the similar magnitude, 
about -0.4 ^B, independently of the alkaline earth ele- 
ment, while that of Re is of the order of -1.1 /is- These 
induced spin magnetic moments are antiparallel aligned 
with the Fe spin moment. Concerning the Fe valency, for 
all the compounds the trivalent configuration has been 
most favourable, followed by the tetravalent and divalent 
ones. 
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FIG. 1: Spin polarized total and oxygen 2p densities of states per formula unit, for Ca2FeRe06 (top row, left), Ca2FeMo06 
(top row, right), Sr2FeMo06 (bottom row, left), and Ba2FeMo06 (bottom row, right). 
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FIG. 2: Spin polarized Fe 3d and Re 5d densities of states per formula unit in Ca2FeRe06 (top row, left), and spin polarized 
Fe 3d and Mo 4d densities of states per formula unit in Ca2FeMo06 (top row, right), Sr2FeMo06 (bottom row, left), and 
Ba2FeMo06 (bottom row, right). 



